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Abstract 

We study Dirac equation in Kerr-Taub-NUT spacetime. We use Boyer- 
Lindquist coordinates and separate the resulting equations into radial 
and angular parts. We get some exact analytical solutions of the an- 
gular equations for some special cases. We also obtain the radial wave 
equations with an effective potential. Finally we discuss the poten- 
tials by plotting them as a function of radial distance in a physically 
acceptable region. 
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1 Introduction 



Kerr-Taub-NUT (Newman- Unti-Tamburino) spacetime is the generalized so- 
lution of the vacuum Einstein field equations in general relativity [1]. The so- 
lution is obtained by introducing an extra non-trivial magnetic mass parame- 
ter which is also called "gravitomagnetic monopole moment" (NUT charge). 
Kerr-Taub-NUT solution describes the spacetime of a localized stationary 
and axially symmetric object [Ij. The solution contains three physical pa- 
rameters: The gravitational mass, which is also called gravitoelectric charge; 
the gravitomagnetic mass that is also identified as the NUT charge; the ro- 
tation parameter that is the angular speed per unit mass. The NUT charge 
produces an asymptotically non-fiat spacetime in contrast to Kerr geome- 
try that is asymptotically fiat [2]. Although the Kerr-Taub-NUT spacetime 
has no curvature singularities, there exist conical singularities on the axis of 
symmetry [3]. One can get rid off conical singularities by taking a period- 
icity condition over the time coordinate. But, this leads to the emergence 
of closed time-like curves in the spacetime. It means that, in contrast to 
Kerr solution interpreted as a regular rotating black hole, the Kerr-Taub- 
NUT solution cannot be identified as a regular black hole solution due to its 
singularity structure. An alternative physical interpretation of Kerr-Taub- 
NUT spacetime can be found in |1] where the NUT metric is interpreted as 
a semi-infinite massless source of angular momentum. Despite the fact that 
Kerr-Taub-NUT solution has some unpleasing properties, it is vastly stud- 
ied for exploring various physical phenomena in general relativity due to its 
asymptotically non-flat spacetime structure O EJ [7] . 

In the present paper, we study the Dirac equation in Kerr-Taub-NUT 
spacetime. Dirac equation has been extensively examined in various grav- 
itational spacetimes including Schwarzschild geometry |8], Kerr-spacetime 
O [TOl [m |12] , 4-dimensional constant-curvature black hole spacetime [13], ro- 
tating Bertotti- Robinson geometry p^, 4-dimensional Nutku helicoid space- 
time [15] and open universe geometry [16] . In all these spacetimes except the 
Kerr geometry, exact solutions of Dirac equation have been found for all pos- 
sible values of physical parameters. However in the rotating Kerr spacetime, 
exact solutions of Dirac equation have been obtained only for some special 
values of the parameters [HI [12] . In [11] , the series solution of angular Dirac 
equation have been given while in [12], angular solutions have been presented 
by using spectral decomposition method in which the angular wave functions 
are expanded in terms of spheroidal harmonics. By this method, a three-term 
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recursion relation is achieved and eigenvalues of the angular equations are 
solved. 

On the other hand, in all these works, the separability of the Dirac equa- 
tion has been accomplished. Separability was first discovered in Hamilton- 
Jacobi and relativistic wave equations by Carter [T7j. Later on it has been 
discussed in higher dimensional spacetimes in the works [ISl HH [20] • Sepa- 
rability also implies the existence of integrals of motion associated with the 
second order symmetric Stackel-Killing tensors fl9[ [20] . On the other hand, 
the separability of the Dirac equation was first discovered by Chandrasekhar 
P, [To] and later on it was extended to Kerr-Newman [21] and AdS-Kerr- 
Taub-NUT spacetimes [2^ . 

In this work, we obtain the set of equations by employing an axially sym- 
metric ansatz for the Dirac spinor. Equations obtained are separated into 
radial and angular parts with appropriate substitutions of spinor fields. We 
try to solve angular equations exactly. But unfortunately, we are unable to 
get exact analytical solutions to general angular equations for all physical 
parameters. Under some restrictions implemented on the separation con- 
stant, we present some exact solutions of the equations with and without 
gravitomagnetic mass and rotation parameters. Indeed, they can be solved 
exactly in terms of hypergeometric functions for the cases where the mass 
of the Dirac particle is equal to or twice the frequency of the spinor wave 
function. In the final part, radial equations are discussed. With some trans- 
formations on the dependent and independent field variables, wave equations 
with an effective potential barrier are obtained. To understand the physi- 
cal behavior of the potentials, they are plotted with changing frequency and 
gravitomagnetic mass parameter in the physically acceptable regions. 

Organization of the paper is as follows: In section 2, we present the 
general form of Dirac equation in exterior forms. In section 3, we obtain Dirac 
equation in Kerr-Taub-NUT spacetime. In the subsections, we discuss the 
separability of the equations, obtain the angular and radial equations. Next, 
we find some exact analytical solutions of the angular equations. Finally, we 
study the radial wave equations and examine the behavior of the potential 
barriers that come out in the transformed radial equations. We end up with 
some comments and conclusions. 
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2 Dirac equation in 4-dimensional spacetime 



We consider a 4-dimensional spacetime manifold M equipped with a Lorentzian 
metric g with signature (—,+,+,+) and a metric compatible connection 
V. We assume that our spacetime manifold has a spin structure group 
5'pm+(3,l). It is known that the fundamental group of Lorentzian group 
SOj^{?>, 1) is 2,2 so that it has a universal covering group of Spin+{3, 1) that 
is the multiplicative subgroup of complex Clifford algebra C/3,1. 
In exterior forms, the Dirac equation can be written as [23] 

*'j ADiIj + ^ip*! = 0, (2.1) 

where 7 is C/s^i-valued 1-form 7 = 7'^ea. We choose the units such that 
c = 1 and h = 1. Here * denotes Hodge-star operator and /i is the mass 
of the particle, {e^j's are the orthonormal co-frame 1-forms such that the 
metric g = r]abe"' ® e^. i/j represents C^-valued Dirac spinor whose covariant 
exterior derivative can be written as 

Di; = dij + ^a'^^'uabiJ, (2.2) 
where a""^ = |[7",7''] and {7"}'s satisfy the relations 

{^^^^} = (^V + ^y) = 2rtx4. 

Uab are the connection 1-forms that satisfy Cartan structure equations 

de" + u%Ae^ = T'' 
where denotes torsion 2-form and metric compatibility implies that Wab = 

-Wba- 

Since is isomorphic to A^4(C) that is the set of 4 x 4 complex 

matrices, we can choose the representation 
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where a* are the Pauli spin matrices and / is the 2x2 identity matrix. 
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3 Dirac Equation in Kerr-Taub-NUT Space- 
time 

In this section, we examine the Dirac equation in Kerr-Taub-Nut spacetime. 
In Boyer-Lindquist coordinates, Kerr-Taub-NUT spacetime can be described 
by the metric with asymptotically non-flat structure, 

g = -^^dt- xdipf + S + de'^ + ^ {adt - (r^ + f + a^)dvY 

(3.1) 

where 

S = + (£ + acos^)^ A = r"^ - 2Mr + - ^^ 

and 

X = a sin^ 9 — 2i cos 9. 

Here, M is a parameter related to physical mass of the gravitational source. 
a is associated with its angular momentum per unit mass and i denotes 
gravitomagnetic monopole moment of the source. For the metric (13.11) . we 
choose the co-frame 1-forms 

e° = f^j (dt-xd^), ""'^Ia) ^^-^^ 



j:'/^d9, = |t7? i^dt - (r^ + f + a^)dip) . 



We consider that the spacetime is Levi-Civita (torsion-free) such that con- 
nection 1-forms a;" f, can be determined from the equation 

de" + uj\ A e'' = 

which has a unique solution 

a;% = ^ {e%L'' tbde,) + n.de" - L^deb) . (3.3) 

Here La = iXa inner-product operators that satisfy L^e'^ = 6^. From 
equation (13. 3p . we can determine connection 1-forms u"" b'- 

w° 1 = A2dt + Asdf, 2 = B^dip, 

ujO 3 = Cidr + C2d9, 2 = Eidr + E2d9, (3.4) 
3 = Gsdip, 3 = K2dt + K^dif, 



where 



A, 



Mr^ + 2e^r + 2lar cos Q - M{i + a cos i 



A, 



X {(m - r)S - 2rr - 2Mr^} - 2ir cos 9{r^ + f + a 

S2 



A^/'^sinOie + acos9) ^ ar sinO ^ A^/^,, 
B3 = ^ Ci = ^^, C2 = ^(£ + acos( 



asm9{i + acosO) rA^/^ rsiiiM^/^ 
= ^175^ , = G, = , (3.5) 



ir^ -{i + a cos 6) (2Mr + f + ai cos 1 
S2 



_ cos 9{r^ + e^ + a^){i: + 2£{e + a cos 9)) + x{2Mr + 2e) (£ + a cos 9) 
Since the space-time is axially symmetric, we can take 



^3(r,^) 



(3.6) 



where m denotes azimuthal quantum number. Then we substitute fl3.4p . 
(13.51) and (13.61) into Dirac equation (12.11) and obtain the following equations: 

[ AV^sing ) SV^ SV^ ^+2(^2+54)^4+2 (51-^3)^4+/.^! = 0, 

(3.7) 

A Vising J ^+EV^ ^ + 2(^4-52)^3+2(5i+53)^3+/.V^2 = 0, 

(3.8) 
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\ A Vising J SV^ ^+SV^ ^+2(^2-54)V'2+2(<^i+<^3)^2+/.^3 = 0, 

(3.9) 

(3.10) 

where 

_ _ (^2 ^ ^2 ^ ^2) ^-^Q ^ ^ 5, _ ^1/2 

_ x^^^'^ + (^^ + + ^^) sin 61 _ a sin 6' + A^/^ 

«3 ^17^ , «4 , 

and 

{I + a cos 9) e ar sin 6 

^1 = ^372 ' 2 = j.3/2 ' 

_ cos ei: - a sin^ 9{l + a cos 6) _ (S(r-M)+rA) 

' ~ sineS3/2 ' - A1/2E3/2 • 

3.1 Separability of the equations 

Adding and subtracting equations < \3.7h . (13. Sp . (13. 9p . (I3.10p and defining 

Fi = + ^2 

^2 = ^2-^1 (3.11) 

^3 = + ^4 
F4 = -04 - ?/;3, 



we obtain 

1 ( ma - uir"^ + + f) z A^/^ 



+ ^ 1 • / + - ^(^ + acos^)) - £ ^ F4 + /iFi = 0, (3.12) 



EV2 1^ AV2 2 S 

1 r m — i a sir 

EV2 I sin0 ^ 2^ 

-ma + a;(r^ + + ^^) i A^/^ 



r-2(^ + acos^)) -V\ F3 



[r -i{i + acose)) - V} 



Si/2 ]^ AV2 2 S 

1 J-m + wx , ^ asin^ . n/Q1Q^ 

SV2 I sing 2^^^^ ~ + acos^)) - £ [> F3 - /xFa = 0,(3.13) 
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^ ^ ^ +-— ^(r + «(£ + acos^))+I?^Fi 



EV2 [ AV2 2 S 

1 f — m + i a sin 6', ^ ^ 

+ ^ 1 sing + 2^^^"" + + «cosg)) + £ [. F2 + ^Fs = 0,(3.14) 



1 f -ma + a;(r2 + + i A^/^ ^, ^ 



SV2 [ Ai/2 2 S 

+^ [-—Q^ + 2^r + + + ^ j i^i - /"i^4 = 0, (3.15) 

where 

At this stage we redefine 

= ^(r-i(£ + acosg))^/^Fl, 

F2 = -i{r -i{e + a cos 9) F2, (3.16) 

F3 = (r + ^(£ + acosg))^/^F3, 

= {r + + a cos 6))^^^ F^. 

Then we multiply equations fl3A2|) and f l333|) by S^/^ _ ^(^^ + acos^))^/^ 
and equations f l3.14p and f l3.15p by S^/^ {r + i{£ + a cos 6))^^"^ and we simplify 
resulting equations. We finally get 

ma - u{r^ + + a^) ^ f rn - ^\ 



Ai/2 j I sin6' 

-li{ir-{i + acos9))Fi = 0, (3.17) 



-m + ux /A ^ f -ma + uj{r^ + e'^ + a^) ( 



sm9 j I A^/^ 

-n{ir-{i + acos9))F2 = 0, (3.18) 

ma-u{r' + e + a') ^ f m - ux A p^ 



Ai/2 j I sing 

+fi{ir + {i + acos9))F3 = 0, (3.19) 
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m — ux 
sin 9 



ma — u{r'^ + £^ + a^) 
A172 

-yLi(«r + (£ + acos6'))F4 = 0. 



(3.20) 



Equations (I3.17p . (13.181) . (I3.19P and (13.201) imply the separability ansatz 

F, = R,{r)S,{e), 

F2 = R2{r)S2{e), (3.21) 
F3 = i?2(r)^i(^), 
F, = R,{r)S2{e). 



With the ansatz above, equations take the following forms: 



ma — ujia^ + + £^1 



Ai/2 



V\R2{r) -ifirRi{r) 



Si{0) 



+ 



m — ux 
sin 6 



C}S2{e) + fi{i + acos9)Si{e) 



Riir) = 0, 



(3.22) 



+ 



—m + ujx 



sin 6 

-ma + a;(a^ + + ^^) 
AV2 



ma — a;(a^ + + £^) 



£^>Si(^)+/i(£ + acos6')52(6') 



i?2(r) 



P ^i?i(r) -i^rR2{r) 



S2{e) = 0, (3.23) 



Ai/2 



I? ^ Ri{r) +ifirR2ir) 



m — tux 
sin ^ 



C} 82(6) + fi{i + acose)Siie) 



R2{r) = 0, (3.24) 



+ 



m — ux 



sin 6 

ma — u{a^ + + ^'^) 
AV2 



+ £ ^ ^1 (^) - ;u(^ + a cos 6')52(6') 



i?i(r) 



P \ R2{r) - inrRi{r) 



S2{e) = 0. (3.25) 
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These equations further imply that 
Aii?i(r) 
A2i?2(r) = 



A3i?2(r) 
\4Ri(r) 





-uia^ + r^ + e) 


1 ma 


Ai/2 


—ma 


+ uia"^ + r"^ + f) 








- uj{a^ + r^ + e) 


1 ma 


Ai/2 


( ma 


-a;(a2 + r2 + £2) 



AV2 



-1 




— ifirRi{r), 


(3.26) 




.i?i(r) 


- iiJ,rR2{r), 


(3.27) 




[i?i(r) 


+ ifirR2{r), 


(3.28) 


+ p 




) - ifirRi{r), 


(3.29) 



and 



X^Si{9) = { -II^,±^ + C} S2{9) - fi{e + acos9)Si{9), (3.30) 
sm fc* 

A2^2(g) = { ""7"^/ +C}Sii9)-fi{i + acos9)S2i9), (3.31) 
sm u 

^3Si{9) = { +c}S2{9)- + a cos 9) Si{9), (3.32) 

sm 9 

X4S2{9) = { ^"^.^ - C } Si{9) + + a COS 9) S2{9). (3.33) 
' sm9 

For consistency of the equations fl3.26p - fl3.33p . we choose Ai = A3 = A4 = A 
and A2 = —A. Then we obtain following independent radial and angular 
equations: 

Aft(.) = + + _ p| _ (3 3„ 

ma — uia^ + + 



and 



Ai?2(r) = { ^jy^ +V\R^{r) + tfirR2{r), (3.35) 



AS'2(6') = j '^X _c}Si{9) + fiii + acos9)S2i9), (3.36) 
sm 

- 1 
sin ^ 



A5i(^)= J '[^ + C\ 32(9)- + a cos 9) Si{9). (3.37) 



3.2 Angular equations 

Angular equations (I3.36P and fl3.37p can be arranged as 
^^^+1(1 + 2w£) cot e-au sin 9+ \ Si{e) = {fi{e + a cos 6) - X) 82(9), 



d9 {'2 ' sin^ 

(3.38) 

and 

^'^2^/(1 _2a;^) cot^ + awsin^ - ^j^sl^) = + acos^) + A) S^{Q). 



d9 y2 ' sin^ 

(3.39) 

By affecting the transformation 

S\ = cos ^0 Ti + sin T2 

82 = - sin ( ^ ) Ti + cos ( ^ ) T2, (3.40) 



2 / V 2 



we get the following equations: 



+ m J cot 9 + i{fi — 2u) sin 9 H — — - + {afi — au) sin 9 cos 9\Ti 



d9 {\2 J / gj^^ 

+ I (m + A + i) + £(2cj - /x) cos 9-au sin^ - cos^ ^ J> = 0, (3.41) 



+ I f ""T- ) cot 9 — — 2uj) sin ^ ^ + (afi — au) sin 9 cos 6* 1 T2 

at/ 1^ \ 2 / sin 6* J 

+ I (m - A - i) + £(2w - //) cos 9-au sin^ - a/i cos^ ^| Ti = 0, (3.42) 

If we redefine Ti = T+ and T2 = T_, we see that f l3.4ip and fl3.42p satisfy 
the following second order differential equations: 

" ^± + + iV±T± = 0, (3.43) 



d9^ " d9 
where 

2 (a// — aa;) sin 6* cos 9 — (2a;£ — /x£) sin 6* 



M± = cot 9 



(=F(| + A) — m) — 1^(20; — //) cos ^ + aa; sin^ 9 + afj, cos^ 6* 

(3.44) 
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and 



[m±l) +Acoseu£ (m±|) + AuH'^ 



+ {±i{fi - 2lo) + (2u;^ - i?)2la) cos 61 + (±2 - a^i - auj){a^ - auj) cos^ 6 
+ {W - T - oj)a - a^u? + 2maw + A(A + 1) (3.45) 

/ (| ± m) cos ^ ± 2a;£ ± sin^ Q - 2u) + - w) cos 0) \ 
\ ("F(i + A) — m) — £(2a; — //) cos 9 + au sin^ 9 + afi cos^ 6' / 

(2a(/i — a;) cos^ — (2a; — . 



Taking x = cos 6*, equations fl3.43p transform into following form: 

doc 

where 



(1 - 1^-)^ + M±^^ + ]V±r± = (3.46) 



(=F (2 + ~ "^j ~ ^(2c<; — + au! + [afi — au)x^ 

and 

(m ± 1)^ + 4a;^ (m ± |) x + 4u;2^2 



A^4 



1 



+ (±£(/i - 2a;) + (2a;2 - ^^)2£a) x + (±2 - a// - aa;)(a/x - 
+ (4cc;2 - ^2)^^ T (^ - oj)a - a^u? + 2macc; + A(A + 1) (3.48) 

/ (i ± m)x ± 2a;£ ± (£(/i - 2a;) + a(/i - u)x) (1 - x^) ' 
\ ("F(i + A) — m) — £(2a; — \i)x + ao; + (a// — au)x'^ 

(2a(/i — a;)x — (2a; — . 
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Now we investigate exact solutions to equations f l3.46p for some special cases: 
i. £ = 0, a = : 

In that case, equations f l3.46p take the simple form, 



A((l-.=)^Zk) + U + i)_(!p^lT,=0, (3,49) 

ax \ ax I \ 1 — x^ 



2' 



In general, the solutions to those equations are associated Legendre functions 
and Q'y^ of the first and second kind respectively which can be expressed 

as 



and 

Ql^ = A\-—-] F( -A,A + l,z/± + l 



x + lj V ' 2 

X + \ \ ' I ^ , I — X 



+B[j^j F(^-A,A + l,l-z/±,^J (3.51) 
with i^-t = m ± I an integer, 

^ ' 2r(l + A-z/±) ' 2 ^ ^ 

When A and i'± are integers {u± being even), solutions describe associated 

Legendre polynomials. 

ii. £ = 0, a 7^ and ^ = u: 

In that case, equations 13.461 can be simplified as 

^((-^=)^)-(^<-^)-^)--«. 

where 

A(A + 1) = A(A + 1) + 2maa; - aV. (3.54) 

In that case, the solutions are the same as in the case (i) except that A is 
replaced by A. 



12 



iii. ^ = 0, a 7^ 0, // 7^ a;: 

In that case, equations 13.461 can be written in the following form: 

(1 _ a;2)^ + (_2x 2iafi-au)x{l-x') \ dT± 

dx"^ I (=F(i + A) - m) + ao; + (a/x - acc;)x2 y dx 

f {m±\Y X 2 / \ 2 2 

+ — h (±2 — a/i — aujjiajj, — au)x =F (/^ — ^jo- — a u 

\ \ — x^ 

+2maa; + A(A + 1) + , ,7 ^ '-^ r± = 0.(3.55) 

(=F(,f + A) — mj + aw + (a/x — aui)x^ 

To simplify the equations I3.55[ we take the constraints 

A = ±afi =1= m — - (3.56) 

in the equations satisfied by r+ and T_ respectively. With these restrictions, 
equations (13.55^ reduce to 

(1 - x' f^ + {Cx' + D^x' + E±) r± = (3.57) 

where 

2 2 2 2 

G = a — a 00 , 

0000 \0 / \ 

D-t = —a/i +2aa; =F'^=Fai^ — — 2aa;m — -, (3.58) 

. ll aO 00"^ 

= — |m±-) ±aa; + A+ 2mauj — a uj — -. 

Equations of the type (13.571) have exact analytical solutions for C = 0, D± = 
or C = 0, D± 7^ [21]. However in our case since fx ^ u, C ^ 0. In 
fact, equations (I3.57P have simple analytical solutions when A = ±fia (which 
corresponds to m = =i=|) which are given by 

T±{x) = d± sin (VCx) + e± cos{VCx) (3.59) 

for n > u, and 

T±(x) = (i±e(^^) + e±e(^^) (3.60) 

for fi < u. 
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iv. £ 7^ 0, a = and /i = 2uj: 

In that case, equations 13.461 take the simple form 

(1 x^)^^^^ 2x^^^ 
dx"^ dx 



1) - ^ -^r^ T± = 0. (3.61) 



+ |A(AH-1)- 



The equation 13.611 can be arranged as 



r/^T dT 
(x^ - 1)'^ + 2x(x=^ - 1)^ + (C± + D±x - A(A + l)x^) T± = 0,(3.62) 

where 

C± = A(A + 1) - (^m± -4a;^£^ D± = -4a;£ (^m ± . 

Under the transformation 

i = \il~ x), y± = (x + l)-^'±(l - x)-'?±T± (3.63) 

the equation satisfied by Y± takes the form 

d^Y dV 
^(^ - + + /5± + 1) - 7±]^ + «± /3± 1± = (3.64) 

whose solution is given by hypergeometric function 

^±(0 = i^(«±,/3±,7±;0- (3.65) 

Hence the solutions are 

T±(x) = (1 + xY^{l - xy^F (^a±, /3±,7±; ^(1 - x)^ . (3.66) 

Here exponents p± and q± are obtained by solving the equations, 

4g±(g± - 1) + 4g± - A(A + 1) + D± + C± = 0, (3.67) 
(P± - g±)(2(p± + g± - 1) + 2) = D±, (3.68) 
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whose solutions are given by 

ql = 1 lm±^ + 2uji] , (3.69) 



1/1 



1 / , 1 



4 V 2 

Pl = 2[rn±^-2u;ej . (3.70) 

a± , f3± and j± can be obtained from 

a± + /3± + l = 2(p± + g± + l), 

«±^± = (P± + q±? -p±-q± + 2(p± + q±) - A(A + 1),(3.71) 
7± = 2g± + l. 

V. £ 7^ 0, a = and fi ^ 2uj: 

In that case, equations 13.461 reduce to 



dx^ \ (:f(| + A) - m) - £(2a; - fi)x I dx 

, (m ± i)2 + 4a;£ (m ± i) X + 4a;2£2 . . . . . 

+ |-^ ±%-2a;)x + A(A + l) (3.72) 

.9 9^ 9 (i±m)x±2w£±£(^-2w)(l-a;2)^ ^ . 

(=F(f + A) - mj - £(2w - /x)x 

As we have done in the case iii., we simphfy f l3.72p by taking the constraints 

fii - 2a;£ = Q + A^ ± m (3.73) 

in the equations satisfied by and T_ respectively. With these restrictions, 
equations (13.721) satisfied for T+ and T_ reduce to 

d'^T dT 
(1 - x^—^ + (1 - x)(l - x')^ + (C+x^ + D+x + T+ = 0, (3.74) 
ax^ dx ' 

and 

d^T dT - 

(1 - x^f—^ - (1 + 2;)(1 - x^)—- + iC.x^ + D^x + E\ T_ = 0, (3.75) 
ax^ ax 
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respectively, where 



D± = -4:ui{m ± 1) - (^m ± , (3.76) 

E± = 4:ui{ui - /i£ ± m) - Q ± {2^1 + 1) . 

Equations of the typ e 13 . 7'^] and [3775] seem harder to obtain for exact analytical 
solutions. However, with X = fii both equations take the following simple 
forms: 

and 

whose solutions are given by 

T+{x) = ci(l + + C2(l + (3.79) 

T_(s) = rfi(l - 2:)('"-i) + (12(1 - x)-("^-5), (3.80) 
where ci, C2, and (i2 are real constants. 



vi. £ ^ 0, a ^ and /i = 2a;: 

In that case, equations 13.461 take the following form: 

d^T± ( 2aujx{l - x2) \ dT-_ 



(l-x')—^-\2x + 



± 



dx'^ \ (=i=(i + A) — m) + aw + awx^ / dx 

+ l-i!!l±i)!±i^i!ll±i)i±i^ + A(A + i) (3.81) 

—Maoj'^x + (±2 — 3auj)aujx'^ =F aw — a^u'^ + 2mau 

, i ± m)x ± 2ui ± aa;x(l — x'^) , 
+2awx \^ ^ ^ ^ 1 T± = 0. 



(t(| + A) - m) + aw + 



atux"^ 
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Under the constraints 

2auj = m± (^ + X] , (3.82) 
equations 13.811 can be simplified as 



(fT 

(1 — x'^Y ^ + (ao + ttix + a2X^ + OaX^ + a^x^) = 0, (3.83) 



and 



1 — X ) ^ + [do + dix + d2X + dsx + d^x j T_ = 0, (3.84) 



dx"^ 

respectively, where 



and 



ai = —Aluo {m+- + auo 

02 = —Qa^u^ + + -"j (^aoj — m — ^ ) , (3.85) 

as = AaluP' ^ 

o 2 2 

a4 = 6a u , 



do = Sa'^u^ - Au^f + ( m - - ) (1 - 2au) 



2 



3 

di = —Aiu { m + auj — - 



a2 = — 6a^aj^ + — ^2aa; — m + , (3.86) 



2 



as = Aaicu 

o 2 2 

a4 = Sa u . 



In order to obtain an exact analytical solution, equation 13.831 can further be 

■I) and £ = ^ 



simplified when A = 2aa; + 1 (which corresponds to m = — |) and " — 



d^T 

[l + x)^^ + (3a^a;2 + 2aa; + 'ia^u'^x) T+ = 0. (3.87) 
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Similarly 13.841 take a simpler form when A = 2au — 1 (which corresponds to 



m = |) and £ = ^ as 



(1 - x)-— + {3a^u^ - 2au - Sa^w^x) T_ = 0. (3.88) 

Both equations are of the type 

dP'v dv 
ic2X + b2)-7^ + (cix + bi)-^ + (coX + bo)y = 0, (3.89) 

with Ci = 0, &i = 0. Under the transformation [2l] 

y = e^-n{z), z = ^{x-Jx) (3.90) 

equation satisfied by VL{z) takes the form 

z^ + ib-z)—-an = 0, (3.91 
az"^ az 

where 

a=— — — , 6 = 0. (3.92) 

2c2k 

Here = ^iid A = — ^. A; can be calculated from C2/c^ + Cq = 0. A 
particular solution of 13.911 (with 6 = 0) can be written in terms of confluent 
hypergeometric function as 

n{z) = z^{a + 1,2, z). (3.93) 

In our case, for T+ 

C2 = 1, 62 = 1, Co = Sa^u'^, 60 = 3a^a;^ + 2auj. (3.94) 
So the solution for T+ can be written as 

T+{x) = e'^'^'^'^z^ia + 1, 2, z) (3.95) 
with z = —2i\/3uja{x + 1) and d = On the other hand, for T_ 

C2 = —1, 62 = 1, Co = —Sa^co'^, bo = 3a^a;^ — 2auj. (3.96) 

In that case the solution for T_ can be written as 

T_(a;) = e-*^"'"^2$(a + l,2,2) (3.97) 

with z = —2i\/3ua{x — 1) and d = —-75- For both and r_, the solutions 
describe oscillating wave solutions with non-uniform amplitudes. 
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3.3 Radial equations 

Radial equations fl3.34p . f l3.35p can be rearranged as 

dRi fr-M .(ma-w(a2 + r2 + £2))2\ 2 , 

^ + [-^ — ^ A J = ^ ^ ^ " 

(3.98) 

dR2 fr-M .{ma-u{a^ + r^ + f))2\^ 2 , „ 

-^+[-^ + ^ ^ ) R. = -^{,r-^X)R,. 

(3.99) 

It can be seen that Ri and R2 satisfy second order differential equations 

^,P,^^,,,^,.0. (3.100, 

where we set Ri = R+ and R2 = R-. Here 

/3(r-M) /^V \^.f M ^ 

= [-^ 7^^' ) I^IvTp ) ' ^'-'^'^ 

and 

L± = Re(L±) +2lm(L±), (3.102) 



where 



Re(L±) = - + — (ma- a;(a2 + r2 + £2)X- _ . 2^2^ (g^^gg) 
(r - M)^V 2/iA 



A(^V2 + A2) A(/iV2 + A 



+ A^,2j2 , ^2^ K« - ^(«' + + ^')) 



and 



1 / r-M \ 

7^I^^-4a;r . (3.104) 



A V/i2r2 + A2 

To get the radial equations in the form of a wave equation, we follow the 
method applied in Chandrasekhar's book [10] but with i ^ 0. Hence we 
consider the transformations 

Pi = lA^/^Ri, P2 = Ai/2^2. (3.105) 
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Let = + + £^ — With these transformations, equations fl3.98p and 
(Kmh take the form 

^^,«i^P..zHMp,. ,3.06, 

(3.107) 

Let 

^ = ^, (3^ = M^-a^ + e. (3.108) 
dr A 

Then in terms of the new independent variable we obtain 

dPi —2i(nr + iX) , /r, , , 

— ^ + 2icoPi = ^ A'/'P2, 3.109 



(3.110) 



where 



2Mr+ + 2£2 - ITH 
w = r H — In 



r — r I 



2/? V ^ 



2^ ln(^^— J, (r>r+) (3.111) 

for > - £2. Here r+ = M + /3 and r_ = M - /3. The relation f l3.11ip 
is single- valued for r > r+ if the inequality 

r\^a^^i^- — >^ (3.112) 

is satisfied. This requires that in the frequency range 

UJ2 < < (x^i, (3.113) 

where 

"^ct ma ^ /o11.^ 

^^= 2Mr^ + 2f2 ' '"^ = FT^' ^'-'''^ 

u becomes single- valued in r in the region where r > r_(_. It is seen that as 
r — 7- oo, w — )■ oo and as r — )■ (r+)+, it — )■ — oo. For completeness, we also 
note that for the critical case when = o? — P' ^ u becomes 

(2£2 + 2M2-i^H) . . . 

M = r-^ ^ + 2Mlnr-M, r>M. 3.115 

r — M 
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'd = arctan (^) . (3.116) 



Concentrating on the case > — i"^, let us make another transformation 

T 

With the new definitions 

Pi = 0le-5*^ P2 = 02e^'^ (3.117) 
equations fl3.109p . f l3.110p take the forms 



Redefining the independent variable as 



4a; V A 

equations (13.118^ and (13.119^ can be simplified as 



u = u arctan , (3.120) 



du 
d(p2 



- + 2iu<pi = W(t)2, (3.121) 
2iujcj)2 = W(pi, (3.122) 



dii 
where 

W= ' , ^ . (3.123) 

(A2 + /.2^2) ^ _ ^ 

By further defining Zi = (pi + 02 and Z2 = 0i — 02, equations f l3.12ip and 
(I3.122P can be rewritten 

dZi 



du 
dZ2 



- WZi = -2iuZ2, (3.124) 

+ WZ2 = -2iuZi. (3.125) 

From equations (13 . 1 24p and fl3.125p we obtain one-dimensional wave equa- 
tions 



+ Au^Zi = V+Zi, (3.126) 



dv? 

dP Zo o 
^ I A, J 



+ Auj'Z2 = V.Z2, (3.127) 



dijL 
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where the effective potentials 



V± = W^± 



dW 



du 



We calculate the potentials as 



P 



where 



2/i2^]r + 2 (A^ + /iV^) r - 

A^A 



2u 



/ = (A2 + ^\^) Q 



(3.128) 



2(a2 + ^V)'/'a1/2 

± 3/iVA ± (A^ + /iV) (r - M) (3.129) 
A^(r - M)' 



(3.130) 



We see that, the effective potentials depend on gravitomagnetic monopole 
moment £ via the functions A and Q, where £ = case is discussed in [TO] . 
We also report that, for yU = 0, the potentials take the simple form 

14(r) = |2AAV2 ± (, _ m) t ^} . (3.131) 



4 Discussion 

Potentials become singular when 7 = 0. Moreover they possess local ex- 
tremums when = which can also be seen in the following graphs. To 
understand the behaviour of the potentials V± in the physical region r > r_|_, 
we plot them as a function of radial distance r for massive and massless 
particles. We take the physical parameters M = 1, A = 1, m = 0.5, a = 0.95. 

Here, Figure [1] and [2] describe the effective potentials V± for massive 
particle with jj, = 1. In the first graph, frequency u is kept constant (w = 0.2). 
We examine the potentials by changing gravitomagnetic monopole moment £. 
We see that, for sufficiently small values of £ including £ = 0, potentials have 
a local maximum in the physical region r > r^. It is also observed that while 
£ increases, local maxima of the potentials seems to disappear. Potentials 
become bounded regardless of the value of £ and approach a constant value 
in the sufficiently large values of r (or r — )■ oo). In the second graph, we 
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12 r 



10 - 




Figure 2: Graph of V± with M = 1, A = 1, /i = 1, ^ = 0.4, m = 0.5, 
a = 0.95 
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Figure 3: Graph of V± with M = 1, A = 1, /i = 0, £ = 0.4, m = 0.5, 
a = 0.95 
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keep gravitomagnetic monopole moment i fixed {£ = 0.4). In that case, we 
investigate the behaviour of the potentials by changing u that can take values 
in the range (13.1131) for -u or m to be single- valued. We see that potentials have 
some local maxima in low frequencies. While the frequency increases, local 
extremums again disappear as in Figure [1] and potentials behave similarly in 
the sufficiently large distances. 

In Figures [3] and HI we plot the potentials for massless (/i = 0) spin- 
I particles. We take same values for the physical parameters. As in the 
previous case, in one plot u is fixed while £ is changed, in the other £ is kept 
constant while u is changed. It can be seen that, contrary to Figures [1] and [2] 
plotted for massive particle, potentials have local maxima for all values of u 
and £. As expected, maximum value of the potentials in the physical range 
r > r+ decreases when compared to potentials in Figures [1] and [21 

5 Conclusion 

In this work, we examine the Dirac equation in 4-dimensional Kerr-Taub- 
NUT spacetime described by the physical parameters, the mass M, angular 
speed a per unit mass and gravitomagnetic monopole moment £. In the ab- 
sence of gravitomagnetic monopole moment £, equations are separable in the 
Boyer-Lindquist coordinates as shown in [21 HO]- We show that, in the same 
coordinate system, equations are still separable in the presence of gravito- 
magnetic monopole moment £. We obtain angular and radial equations for 
arbitrary £. We find some exact solutions to the angular equations with and 
without gravitomagnetic monopole moment £ and rotation parameter a. We 
see that when the mass of the particle is equal to or twice the frequency of 
the spinor wave function, some angular solutions can be represented in terms 
of hypergeometric functions. 

We also discuss the radial equations and get a wave equation with an 
effective potential. To realize the physical interpretations of the potentials, 
we plot them as a function of radial distance r. From the plots, it can be 
seen that the strength of the potential barrier decreases while NUT charge 
£ increases. We believe that to better understand the physical significance 
of the NUT charge, Dirac Hamiltonian should be constructed for the Dirac 
equation [23]. This will be the subject of our future research. 
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